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ABSTRACT 

We have investigated the mass-metallicity (M-Z) relation using galaxies at 0.4 < z < 1.0 from the 
Gemini Deep Deep Survey (GDDS) and Canada- France Redshift Survey (CFRS). Deep K and z' band 
photometry allowed us to measure stellar masses for 69 galaxies. From a subsample of 56 galaxies, 
for which metallicity of the interstellar medium is also measured, we identified a strong correlation 
between mass and metallicity, for the first time in the distant Universe. This was possible because of 
the larger base line spanned by the sample in terms of metallicity (a factor of 7) and mass (a factor 
of 400) than in previous works. This correlation is much stronger and tighter than the luminosity- 
metallicity, confirming that stellar mass is a more meaningful physical parameter than luminosity. 
We find clear evidence for temporal evolution in the M-Z relation in the sense that at a given mass, 
a galaxy at z ~ 0.7 tends to have lower metallicity than a local galaxy of similar mass. We use the 
z ~ 0.1 Sloan Digital Sky Survey M-Z relation, and a small sample of z ~ 2.3 Lyman break galaxies 
with known mass and metallicity, to propose an empirical redshift-dependent M-Z relation, according 
to which the stellar mass and metallicity in small galaxies evolve for a longer time than in massive 
galaxies. This relation predicts that the generally metal poor damped Lyman-a galaxies have stellar 
masses of the order of 10 8 8 M Q (with a dispersion of 0.7 dex) all the way from z ~ 0.2 to z ~ 4. 
The observed redshift evolution of the M-Z relation can be reproduced remarkably well by a simple 
closed-box model where the key assumption is an e-folding time for star formation which is higher 
or, in other words, a period of star formation that lasts longer in less massive galaxies than in more 
massive galaxies. Such a picture supports the downsizing scenario for galaxy formation. 
Subject headings: galaxies: fundamental parameters - galaxies: abundances - galaxies: ISM - ISM: H 
II regions - galaxies: evolution - cosmology: observations 



1. INTRODUCTION 

Our exploration of the evolution of the cosmic metal 
enrichment relies mainly on two methods. One is based 
on the detection of absorption lines in the neutral inter- 
stellar medium (ISM) of galaxies crossing QSO sight-lines 
(Prochaskaet al. 2003), and gives information of one line 
of sight in the galaxy. The other uses emission lines of 
the warm ISM (HII regions) detected in the integrated 
galaxy spectra, and is observationally much more chal- 
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lenging at high redshift because important lines get very 
weak and redshifted to the near IR. 12 

The warm ISM metallicity and the mass of galaxies 
are strongly correlated in the low-z Universe (Lequeux 
et al. 1979). More massive galaxies have higher metal- 
licity than less massive galaxies. This mass-metallicity 
(M-Z) relation has never been detected at high redshift. 
In fact, measurements of the stellar mass (strongly cor- 
related with the dynamical mass; Brinchmann & Ellis 
2000) require deep optical/NIR photometry of faint tar- 
gets, and are not easy to obtain for a sufficiently large 
sample for which metallicity is known. 

To explore any evolution with redshift, luminosity, as 
a proxy for mass, can be used instead, although lumi- 
nosity is more difficult to interpret physically than mass. 
The metallicity of large 0.3 < z < 1 galaxy samples has 
been studied in relation to the galaxy luminosity by Kob- 
ulnicky et al. (2003), Lilly, Carollo, & Stockton (2003), 
Kobulnicky & Kewley (2004) and Liang et al. (2004). 
Kobulnicky & Kewley (2004; hereafter KK04) have used 
about 200 galaxies from the Team Keck Redshift Sur- 
vey (TKRS; Wirth et al. 2004) and found a luminosity- 
metallicity (L-Z) relation displaced towards higher lu- 
minosities with respect to the same relation detected at 
z = (Kennicutt 1992; Jansen et al. 2000). This sug- 
gests a redshift evolution of the L-Z relation. At higher 
z, only a handful of galaxies have been studied so far 

12 Stellar features can also be used as metallicity estimator, how- 
ever this way it is much harder to obtain accurate information for 
a single galaxy (Pettini et al. 2000; de Mello et al. 2004). 
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(Pettini et al. 2001; Shapley et al. 2004) and again low 
metallicities are found compared to galaxies with similar 
luminosity at low redshifts. 

Do we also expect the M-Z relation to evolve with 
time? Shifts of the L-Z relation will arise from sim- 
ple luminosity evolution (i.e. changing the mass-to- 
light ratio) and this is well documented in galaxies over 
< z < 1 (Schade et al. 1996; Vogt et al. 1996). In 
particular passive fading of stellar populations (as ob- 
served in early-type galaxies, e.g. Aragon-Salamanca et 
al. 1993) would shift L-Z but not M-Z. However, astro- 
physically, we would expect the star-formation activity 
over < z < 1 to manifest itself in changes of metallicity 
as well as mass. The question is what is the form of this 
change and can it be detected? 

From a different prospective, the cold ISM, damped 
Lyman-a (DLA) galaxies detected in QSO spectra have 
metallicities ~ 1/10 solar at z ~ 0.7 (Prochaska et al. 
2003), i.e. one order of magnitudes lower than line emit- 
ters at similar redshifts (KK04). Part of the discrepancy 
could be due to the fact that the absorption lines statis- 
tically probe the outskirt of galaxies, where metallicity 
is lower, while emission line flux originates in the cen- 
tral more metal rich region (Chen, Kennicutt, & Rauch 
2005; Ellison et al. 2005; Schulte-Ladbeck et al. 2005). 
Another possibility is that DLAs and line emitters are 
probing different galaxies all together. Unfortunately, so 
far no information on masses of DLA galaxies is available. 

In this work we present the first attempt to derive 
the mass-metallicity relation in high redshift galaxies se- 
lected by near-IR photometry down to small masses. The 
sample is selected from the 0.4 < z < 1 galaxies of the 
Gemini Deep Deep Survey (GDDS; Abraham et al. 2004) 
and Canada-France Redshift Survey (CFRS; Crampton 
et al. 1995; Le Fevre et al. 1995; Lilly et al. 1995). Masses 
were estimated using deep optical-NIR photometry and 
metallicities using optical spectroscopy. For consistency 
with GDDS, metallicities for the CFRS galaxies were re- 
computed using emission line fluxes published by Lilly et 
al. (2003; hereafter LCS). 

Before continuing our discussion it is important to em- 
phasize that measurements of physical parameters can 
differ by some factor, depending on how they are esti- 
mated. For instance, masses can change by a factor of 
~ 2, depending on which initial mass function (IMF) 
is applied. Similarly, metallicities can be different by a 
factor of 2 — 3, depending on which set of lines and/or 
calibrator are considered (Kennicutt, Bresolin & Garnett 
2003; KK04). Thus, we paid special attention to deriv- 
ing results, and comparing them with other works, in a 
consistent fashion. 

The paper is organized as follows: in §2, we describe 
the sample section from the GDDS; in §3 we present the 
GDDS composite spectrum used to derive an overall stel- 
lar absorption and dust extinction; §4 & §5 describe the 
stellar mass and the metallicity derivations, respectively, 
for the GDDS and the CFRS samples; §6 & §7 are about 
the luminosity-metallicity and mass-metallicity relations, 
respectively; in §8 we discuss possible systematic effects; 
§9 is on the modeling of the redshift evolution of the 
M-Z relation; the discussion and the concluding remarks 
are in §10 and §11. Throughout the paper we adopt a 
h = H o /100 = 0.7, n M = 0.3, Q A = 0.7 cosmology 
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Fig. 1. — The Iab magnitude vs. redshift of the GDDS galaxies 
at 0.40 < z < 0.98. The selected galaxies are shown as filled 
squares; 25% and 75% are classified as late-type and intermediate- 
type galaxies, respectively. For the rest of the sample, pure early- 
type and intermediate-type galaxies are shown as empty circles and 
squares, respectively. The non-selected intermediate-type galaxies, 
generally show only [Oil] emission and no other line. One galaxy 
(SA22-2107) is not included in the selected sample because it is 
contaminated by an interloper. Crosses are galaxies not included 
in the sample because one or more emission lines are corrupted or 
outside the spectral range. 



(Spergel et al. 2003). 

2. THE GDDS SAMPLE SELECTION 

The Gemini Deep Deep Survey 13 (Abraham et al. 2004) 
is the deepest survey targeting galaxies in the redshift 
desert (0.8 < z < 2). It is based on spectra obtained 
with the Gemini Multi Object Spectrograph (CMOS), 
operating in nod & shuffle mode (Glazebrook & Bland- 
Hawthorn 2001; Cuillandre et al. 1994) for precise sky 
subtraction. The survey has been primarily targeting 
galaxies with photometric redshifts z > 0.8 and K < 

20.6, but other objects were included to fill gaps in the 
CMOS masks. The galaxy sample is complete down to K 
(Vega) magnitudes K = 20.6, and I magnitudes Iab — 

24.7. For comparison, the Iab magnitude limit for the 
TKRS of KK04 and CFRS of LCS is I AB » 24 and 22.5, 
respectively. For the z < 1 galaxies, GDDS is much 
deeper than CFRS because it was designed to reach L* 
for local galaxies to z — 2, the same flux limit is sub-L* 
at z < 1. The observed GDDS I magnitude corresponds 
to the rest-frame B magnitude for a galaxy at z ~ 0.7. 

In terms of stellar mass the GDDS is complete (for 
z < 2) down to M* = 10 10 8 M Q for all galaxies (Glaze- 
brook et al. 2004). The completeness increases to 
M* = 10 101 Mq and 10 9 6 M , for star-forming galaxies 
at z < 2 and z < 1, respectively. More than 300 spec- 
tra were taken, of which about 200 are z > 0.3 galaxies 
with secure redshifts. The slit aperture used for the ob- 
servations was 0.75 x 1.1 arcsec 2 , and the signal in the 
direction perpendicular to the dispersion was extracted 
over ~ 0.8 arcsec. The wavelength range observed is 
typically 5500-9800 A. 

For this work, the GDDS sample selection is based 
on the requirement that the spectrum of the galaxies 
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(with secure redshift) covers the spectral interval of the 
[OII]A3727, [OIII]AA4959,5007 and H/3 lines. These are 
the lines necessary to determine the metallicity using the 
R23 calibrator (Pagel et al. 1979). The redshift range is 
thus restricted to 0.40 < z < 0.98. Figure Q shows Iab 
vs. z for GDDS galaxies in this redshift interval (73). 
Twelve of these are excluded from the sample because 
one or more lines are corrupted or just outside the spec- 
tral range (which to some extend depends also on the 
position of the slit in the mask); one object is clearly 
contaminated by an AGN and another by an interloper, 
so are excluded from the sample. Among the 59 remain- 
ing objects, the selected 28 (filled squares) show [Oil], 
H/3 and [OIII] emission; 25% and 75% are classified as 
late type and intermediate type, respectively. The other 
non-selected galaxies are early type (empty circles) or in- 
termediate type (empty squares) with only weak or no 
[Oil] emission, and no other emission line. 

Our detection limit for emission line fluxes is a function 
of redshift (the spectral sensitivity is higher in the blue 
than in the red) and is well represented in the observed 
redshift interval by the function 

JW10- 18 ergs" 1 cm" 2 ) = 4.31-z - 1.12 (1) 

We detected line fluxes down to (0.G-3.2) x 10~ ls erg s _1 
cm~ 2 in the interval z = 0.4 — 1 (3a). For 4 galaxies, the 
[OIII] or H/3 emission is below detectability, so a limit was 
derived for the metallicity. The CFRS detection limit of 
LCS is - 2 x 10" 17 erg s' 1 cm" 2 , i.e. a factor of 6 - 30 
higher than our flux limit. 14 

Spectra of the selected sample are shown in Figures |2- 
[SJ For 15 of the galaxies we also obtained Hubble Space 
Telescope (HST) Advances Camera for Survey (ACS) im- 
ages (Figure EJ). For the majority of them a disk struc- 
ture is apparent. To measure fluxes, the continuum was 
estimated from the mean value in two small spectral re- 
gions before and after the line. Errors are derived using 
the noise spectrum, and generally agree (within a factor 
of 1.5) with the pixel-to-pixel standard deviation in the 
signal spectra. For about 1/3 of the galaxies, the H7 
line is also detected. All line fluxes, not corrected for 
slit aperture losses, are reported in Table 1. The spec- 
tra of SA22-2541 and SA15-4662 show spurious lines due 
to order overlap in the GMOS mask (see Abraham et al. 
2004). Order overlap has negligible effect on the emission 
flux, therefore metallicity, measurements. 

The multi-band photometry for the GDDS is pro- 
vided by the Las Campanas Infrared Survey (LCIRS; 
McCarthy et al. 2001; Chen et al. 2002) and is gener- 
ally accurate (errors are generally less than 1/10 mag; 
Abraham et al. 2004). K and z 1 band magnitudes (Ta- 
ble 1) are key parameters used to measure the stellar 
mass (Glazebrook et al. 2004). The rest- frame -B-band 
absolute magnitude is derived by approximating the SED 
between the V, R and / bands with a power law. 

3. THE GDDS COMPOSITE SPECTRUM 

We created a composite spectrum from all galaxies in 
the GDDS sample, for two reasons. One is to estimate 
the stellar absorption used to correct the Balmer emis- 
sion lines; the other is to estimate an average value of 

14 The flux limit for TKRS is not available, because spectra are 
not flux calibrated. 



dust extinction, via the Balmer decrement. Fluxes of 
Balmer emission lines in Table 1 are only corrected for 
the stellar Balmer absorption, and not for the dust ex- 
tinction. 

3.1. Stellar Balmer absorption 

To combine all spectra together, each spectrum has 
been normalized to unity, by dividing the flux by its 
mean value in the interval AA = 4200 - 4400 A. Thus 
each galaxy has the same weight in the composite spec- 
trum (Figure 0). The best fit to the stellar continuum, 
obtained using Bruzual & Chariot (2003) stellar popula- 
tion synthesis models, is for a 50 Myr old stellar popula- 
tion and a visual extinction 15 for the stellar continuum 
Ay = 1.6. A Calzetti extinction law (Calzetti 2001) and 
a solar metallicity are assumed. The model gives an ac- 
ceptable fit when the dust extinction and the age of the 
stellar population span 0.7 < A v < 2.0 and 30 — 200 
Myr, respectively. If the dust extinction is in the higher 
or lower end, only ages below 80 Myr or in the range 
100 — 200 Myr are acceptable, respectively. In these in- 
tervals, the emission fluxes in the composite generally 
vary by less than 5%. The best-fit model is only used 
for the Balmer absorption correction and does not give a 
comprehensive description of the underlying stellar pop- 
ulation. In general, the properties (e.g., dust extinction, 
metallicity and ages) of the stellar component in the op- 
tical (dominated by small and intermediate mass stars) 
can differ from those of the young star-forming compo- 
nent represented by the optical emission lines. Moreover, 
strictly speaking, the true Ay is very likely lower than 
the estimated value, because of blue-flux spectral loss 
due to the atmospheric dispersion. This has little effect 
on the Balmer absorption features. 

From the composite we derive an equivalent width EW 
correction for H/3 and H7 absorptions of the order of 3.6 
A and 3.4 A, respectively, depending on the underlying 
emission line FWHM. Changing the H/3 absorption EW 
by ±%50 changes the metallicity by less than 0.1 dex in 
83% of the galaxies, and on average by 0.06 dex. For 
the CFRS sample, LCS have used a Balmer absorption 
correction of 3 A, with 2 A uncertainty, as derived for 
local irregulars and spirals by Kobulnicky, Kennicutt, & 
Pizagno (1999). In extragalactic HII regions the correc- 
tion is lower, with a typical EW of 2 A (McCall, Rybski, 
& Shields 1985). 

After subtracting the stellar continuum (bottom panel 
of Figure EJ) , emission line fluxes in the composite are 
measured and reported (relative to H/3) in Table 2. Er- 
rors are estimated by using the standard deviation in the 
subtracted spectrum around each emission line. 

3.2. Balmer decrement and dust extinction correction 

The comparison between the Balmer decrement (high 
order Balmer emission line fluxes relative to the H/3 flux, 
after stellar absorption correction) observed in the com- 
posite spectrum and the theoretical value (in the case 
of no extinction, from atomic physics) provides a rea- 
sonable estimate of the mean dust extinction. The the- 
oretical line ratios (Osterbrock, 1989) are only a weak 

15 Throughout the paper Ay is the visual extinction of the stellar 
continuum, and Ay is the visual extinction affecting the Balmer 
emission lines. 
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Fig. 2. — The sample of GDDS galaxies at 0.4 < z < 1. Relevant emission lines are marked by the vertical dashed lines. 
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Fig. 3. — As in Figure 121 The spurious emission line at A ~ 4000 A in the spectrum of SA15-4662 is due to order overlap in the GMOS 
mask (see text). 
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Fig. 5. — As in Figurc[5] The two spurious emission lines at A ~ 4200 A in the spectrum of SA22-2541 are due to order overlap in the 
GMOS mask (see text). 
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Fig. 6.— HST/ACS F814W images of 15 of the 28 GDDS selected galaxies. The size of the images is 5 X 5 arcsec 2 . 



function of the gas temperature, and are shown in Ta- 
ble 2 for T = 10, 000 K and 5000 K. The visual ex- 
tinction is estimated assuming the MW extinction law 
(practically unchanged for LMC or SMC extinctions). 
The error- weighted mean, from H7/H/? through H8/H/3, 
is A v = 2.13 ± 0.32 or 1.92 ± 0.32 for T = 10, 000 K or 
5000 K, respectively. 

In principle the assumed visual extinction can be 
checked by comparing line EW ratios (almost unaffected 
by dust) with the extinction-corrected flux ratios. This 
method is not very accurate because it is based on the 
assumption that the continuum flux ratio at the [Oil] 
and H/3 wavelengths is nearly one (Kobulnicky & Phillips 
2003). Nonetheless, in the subsample of 15 galaxies for 
which an accurate estimate of the emission line EWs is 
possible, we found a good match between EW ratios and 
flux ratios for Ay ~ 2. 

We compare our Ay with other estimates derived for 
different samples. Cid Fernandes et al. (2005) found 
an empirical relation between the visual gas extinction 
Ay and the visual stellar extinction Ay using data from 
the Sloan Digital Sky Survey (SDSS). If transformed to 
take into account a small difference in the extinction law 
adopted by us, this relation is: 

Ay = 3.173 + \MlAy - 6.418 log f §E ) (2) 

where (Ha/H/3) t h is the line flux ratio expected from 
atomic physics theory. Considering that from the com- 
posite spectrum best fit we derived an upper limit for the 



stellar extinction of Ay = 1.6, then Ay < 3.2 or < 3.0, 
for T = 10, 000 K or 5000 K, respectively We also con- 
sider that Eq.|21is based on data at z ~ 0.1 and can vary 
in the high-z Universe. 

An extinction for the ISM similar to ours (Ay 
2.-1 ■ 0. 1 was derived from the Ha/H/3 value (after stel- 
lar absorption correction) measured in the SDSS com- 
posite spectrum, the "cosmic optical spectrum" (Glaze- 
brook et al. 2003). LCS assumed Ay = 1, and this could 
be more appropriate for relatively brighter (hence likely 
less extincted) galaxies as in the CFRS sample. This 
gives H7/H/3 = 0.39, which is significantly larger than 
our H7/H/3 = 0.29 ± 0.02 (Table 2). An older and more 
extincted stellar population model (older than 100 Myr, 
and Ay > 1.5) would give a higher H7/H/3 in our com- 
posite, then a lower Ay. However, the fit is much worse, 
the x 2 is 3 times higher. 

Calzetti (1997) found a mean visual extinction in a 
sample of 19 starburst galaxies of Ay = 1.35 with a la 
dispersion of 0.77 magnitudes. The higher value in our 
sample can be an indication of a higher dust content in 
galaxies that are generally more massive and more metal 
rich than the local starbursts (Brinchmann et al. 2004) . 

The H7 emission line is detected in 10 GDDS galaxies 
(Table 1), so the Balmer decrement provides the dust ex- 
tinction in these spectra individually. For this subsample 
we found < Ay >= 1.66 with a dispersion of 1.3 magni- 
tudes, indicating that galaxies for which H7 is detected 
tend to be less extincted. However, because the errors are 
always higher than 0.5 mag, we ignore Ay estimated this 
way. It is worth noticing that the median stellar mass 
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for these galaxies (M+ = 10 9 ' 4 M Q ) is lower, but close, 
to the median value for the whole sample (M+ = 10 9 6 
M Q ), suggesting that if the dust extinction were larger 
for more massive galaxies, this effect is not apparent in 
our intermediate redshift sample. 

The visual extinction uncertainty (0.3 mag) does not 
affect very much the metallicity estimates from line ra- 
tios because the wavelength baseline of the emission lines 
used (from [Oil] to [OIII]) is not very large (3727 - 5007 
A). Consequentially our results are relatively unaffected 
for small variations of Ay. We will adopt Ay — 2.1 ±0.3 
and discuss results for Ay = 1. In the GDDS sample the 
metallicity decreases on average by 0.2 dex if Ay goes 
from to 3. 

4. STELLAR MASSES 

The stellar mass of galaxies in the GDDS and CFRS 
samples have been estimated using the procedure de- 
scribed in Glazebrook et al. (2004) and also used in 
Juneau et al. (2005). The galaxy SED is modeled us- 
ing the multi-band photometry. Galaxy model spectra 
are provided by PEGASE.2 (Fioc & Rocca-Volmerange 
1997, 1999). The fitted parameters are the stellar dust 



extinction (in the range < Ay < 2), the age of the stel- 
lar population, and the metallicity. The star formation 
rate is assumed to decline exponentially, with e-folding 
time ranging from r — 0.1 to 500 Gyr, combined with a 
bursty component with r = 100 Myr and variable mass 
fraction. The best-fit SED gives a M/Lk, from which 
masses are reliably derived. The IMF used is that de- 
rived by Baldry & Glazebrook (2003; hereafter BG03), 
and gives a total stellar mass 1.8 or 1.2 times smaller 
than for a Salpeter or a Kroupa IMF, respectively. The 
difference with Salpeter IMF is mainly in the low mass 
stars. The error in the stellar masses includes uncertain- 
ties on Ay, t and mass fraction of the burst. Our mass 
estimates, based on the rest-frame optical-NIR photome- 
try, is much more robust than when based on rest-frame 
UV light, which is strongly affected by very uncertain pa- 
rameters like dust extinction and recent star formation. 
This SED-fitting method gives a mass accuracy gener- 
ally better than a factor of two for the z > 0.8 GDDS 
galaxies (Glazebrook et al. 2004). A similar method was 
recently tested against other stellar mass estimates, and 
was found to be generally robust and accurate (Drory, 
Bender & Hopp 2005). 
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Fig. 8. — Stellar Mass and B absolute magnitude (or luminos- 
ity, upper horizontal axis) for the GDDS (filled circles) and CFRS 
(empty squares) z < 1 galaxies. The open triangles are the z ~ 2.3 
LBGs (Shapley et al. 2004) after correcting the stellar mass for 
BG03 IMF. The straight line marks a constant ratio of stellar mass 
to B luminosity, in solar units, of M^q/Lb q = 0.3. 
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Fig. 9. — -R23 vs. O32 for the GDDS (filled circles) and the 
CFRS (open squares) galaxies. The dust correction is not applied. 
The CFRS points are calculated using emission line fluxes from 
LCS. The dashed line marks the region with O32 = 1. The large 
fraction of galaxies with O32 > 1 in GDDS with respect to CFRS 
is an indication of lower metallicity, and/or higher dust extinction, 
and/or higher ionization. The O32 value decreases by 0.4 dex if 
the dust correction goes from Ay = to 3. 



For the GDDS galaxies in our z < 1 sample, the VIz'K 
photometry allowed reliable mass estimates to be derived 
for 27 galaxies, out of 28 (Table 1). For one galaxy, the 
stellar mass could not be calculated for incomplete photo- 
metric information. The median error is 0.3 dex, and the 
errors are lower than a factor of two and three for 55% 
and 94% of the objects, respectively. For many galax- 
ies, only upper limits to the K magnitude are available, 
therefore stellar masses rely on the z' magnitude. These 
masses were not derived in Glazebrook et al. (2004), but 
their uncertainties are still acceptable, as the observed z' 
band corresponds to the rest-frame V band at redshift 
~ 0.7. These errors are typically ~ 0.15 dex higher than 
those for galaxies with measured K magnitude. In three 
cases for which neither K nor z' are available, mass errors 
are larger. 

The sample median redshift and redshift interval are 
z = 0.79 and 0.47 < z < 0.96, respectively. The lowest 
and the highest masses are M* = 10 8 2 M Q and 10 107 
Mg, and the sample is complete down to M± = 10 96 
M . 

For the CFRS sample of 69 galaxies, VIK photometry 
is used and masses are measured for 42 galaxies (Table 3). 
The mass error is always smaller than 0.6 dex (less than a 
factor of 4), and the median error is 0.3 dex. The lowest 
measured mass is M* = lO 9 4 M Q , i.e. 16 times higher 
than the lowest stellar mass in the GDDS. The highest 
mass is 10 10 ' 8 M©. The median redshift and interval are 
z = 0.68 and 0.479 < z < 0.915, respectively. 

The total number of galaxies for which we measured 
masses is 69. For 56 of these we also determined metal- 
licities (see [0. 

Figure |H1 shows the stellar mass and the B absolute 
magnitude (or luminosity) for the GDDS and CFRS sam- 
ples. It is apparent that GDDS galaxies are less massive 
(on average 2.24 times) and fainter (~ 1.75 magnitudes, 
or a factor of 5 in flux) than CFRS galaxies. GDDS 
galaxies are on average > 1.1 mag fainter in K than 
CFRS galaxies. The fact that the ratio between the mean 
CFRS B luminosity and the mean GDDS B luminosity is 
larger than the mass ratio, indicates that GDDS galaxies 
are on average more dust extincted than CFRS galaxies. 

As a reference, the LBGs at z ~ 2.3 (Shapley et al. 
2004, also in Figure are much brighter and more mas- 
sive than the CFRS galaxies. (LBG masses are scaled 
to take into account the different IMF used, and the re- 
cycling into the ISM, which Shapley et al. did not in- 
clude and which makes the final stellar masses 1.4 times 
smaller.) 

5. THE METALLICITY USING THE R 23 PARAMETER 

Emission line fluxes of [OIIJA3727, H/3, and 
[OIIIJAA4959, 5007 provide an estimate of metallicity 16 
through the R23 calibrator. When the [OIII]A4959 emis- 
sion line is barely detected, the flux is assumed from 
atomic physics to be 0.34 times the [OIIIJA5007 flux, and 
the error propagated accordingly. The R23 parameter, as 
defined by Pagel et al. (1979), is 

, „ , / /[Oil] A3727 + /[OKI] AA4959, 5007 \ /0 n 
l0gi? 23 = fog y~ ■ J = X (3) 

16 We express metallicity in terms of the number of oxygen atoms 
relative to hydrogen: 12 + log(0/H). 
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Fig. 10. — Metallicity as a function of redshift for GDDS (filled 
circles), CFRS (empty squares) and TKRS (crosses, from KK04). 
The mean metallicities in the GDDS, CFRS and TKRS are 12 + 
log(0/H) = 8.63 ± 0.26, 8.78 ± 0.17 and 8.72 ± 0.17, respectively. 
The mean total metallicity (dashed line) is 12 + log(0/H) = 8.72± 
0.17, i.e. roughly solar. We do not detect any significant redshift 
evolution for the total sample of 261 galaxies, from z = 0.32 to 
z = 0.96 (time interval of 4.3 Gyr). 



We also define the O32 parameter: 



log 32 = log 



/[OIII] 



AA4959.5007 



/[OII]A3727 



y 



(4) 



used to correct for the ionization. The R23 and O32 for 
GDDS galaxies are reported in Table 1 and are corrected 
for dust extinction. When dust extinction is applied, R23 
increases and O32 decreases. O32 decreases by 0.4 dex if 
one corrects for Ay in the range — 3. 

The ionization-sensitive O32 not corrected for dust in 
GDDS and CFRS (Figure gj is > 1 in 62% & 27% of 
the galaxies, respectively. In the ~ 200 TKRS galaxies, 
KK04 found that 32% have O32 > 1. However, because 
KK04 used equivalent widths (not very sensitive to dust 
extinction) this fraction would be higher if unextincted 
fluxes are used instead. As noted by Kewley & Dopita 
(2002), 32 is higher not only in highly ionized or dusty 
regions, but also in low metallicity regions. In fact, in 
high metallicity regions, most of the nebular cooling oc- 
curs through far-IR oxygen lines, so the [OIII] emission 
becomes weak and the 032 parameter decreases. 

The generally higher O32 in the GDDS with respect 
to CFRS (Figure EJ) suggests that GDDS galaxies are 
more ionized and/or more extincted and/or less metal 
rich than CFRS galaxies. Additional evidence for more 
extinction in GDDS was provided in § 0] 

The LBGs studied by Pettini et al. (2001) all have 
032 > 1 an d R23 > 5. The difference between these and 
our z r~ 0.7 galaxies likely implies lower metallicities for 
LBGs, not surprising at higher redshifts. 

5.1. The GDDS and CFRS metallicity 

The average metallicity of galaxies can be derived 
through different calibrators which use i?23 and 032- Ini- 
tially proposed by Pagel et al. (1979), the i?23-metallicity 



relation was reformulated by several authors with the aid 
of the ionization-sensitive parameter O32 (Edmunds & 
Pagel 1984; Dopita & Evans 1986; McGaugh 1991; Pilyu- 
gin 2001; Kewley & Dopita 2002; KK04). The difference 
between the different calibrators can be as high as 0.2 
dex, with the theoretical methods giving generally higher 
metallicities than the empirical methods (KK04). An- 
other limit of the R23 calibrator is its double solution. In 
fact, at low metallicity, when the electron temperature T e 
is high enough to keep the gas collisionally ionized, i?23 
scales with metallicity. But when the metallicity reaches 
a limit, the gas cooling occurs through the far-IR oxygen 
lines, then the i?23-metallicity dependence turns around, 
i.e. i?23 decreases for higher metallicity. This happens 
when i?23 is of the order of 10 and 12 + log(0/H) ~ 8.3. 
In the turnover region, the metallicity is easily 0.2 — 0.4 
dex more uncertain because the two solutions are close 
and both possible. To break the degeneracy, other meth- 
ods requiring other lines can be used, for instance [Nil] 
and [SIII] (Denicolo, Terlevich, & Terlevich 2002; Pettini 
& Pagel 2004; Bresolin et al. 2005). However, at high 
redshift this requires expensive NIR spectroscopy, not 
easily accessible. Because T e is sensitive to metallicity, 
the [OIII] A4363/ [OIII] A5007 temperature-sensitive ratio 
gives the most reliable method (Kobulnicky et al. 1999). 
However [OIII]A4363 is a weak line, useful only when the 
spectral S/N is high. 

When the only available set of lines is [Oil] , [OIII] and 
H/3, as in our and all cases where optical spectroscopy 
is available for high-z galaxies, the mass and the age of 
the galaxy can be used to break the degeneracy, as mas- 
sive and evolved galaxies are typically metal rich. As the 
galaxies in our sample are generally massive, and our i? 2 3 
is generally far from the turnover region, we will only 
consider the upper branch solution. This is supported 
by results on galaxies with similar luminosities and red- 
shifts in TKRS and CFRS, for which other metallicity 
diagnostics indicated that the upper branch solution was 
correct (KK04, LCS). 

Among the different calibrators, the empirical one 
proposed by McGaugh (1991) is the most commonly 
adopted. KK04 compared some of them, discussed possi- 
ble biases, and concluded that the best one is represented 
by a mean value between McGaugh's and the theoret- 
ical one given by Kewley & Dopita (2002). The new 
calibrator is approximated, for metallicities higher than 
12 + log(0/H) = 8.4, by the following expression: 

12 + log(0/H) = 9.11 - 0.218x - 0.0587x 2 - 0.330a; 3 - 
+0.199a; 4 - y(0.00235 - 0.1105a; - 0.051a; 2 - 0.04085x 3 - 

+0.003585a; 4 



(5) 



Roughly speaking, the metallicity of Eq. is not too 
different from McGaugh's, on average 0.07 dex higher in 
the range 8.4 < 12 + log(0/H) < 9.3. 

For this work, we used Eq. [S] For 4 galaxies, i?23 is in 
the turnover region, and the metallicity is more uncer- 
tain. In these cases (when 12 + log(0/H) < 8.4) we used 
the McGaugh calibrator and applied the small positive 
correction of 0.07 dex (as derived from the comparison 
with KK04). These galaxies do not change much our 
findings on the M-Z distribution at z ~ 0.7 (see JJ. 
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GDDS metallicities are reported in Table 1 (calculated 
for Ay — 2.1 ±0.3). The mean value (and la dispersion) 
is 12 + log(0/H) = 8.63 ± 0.25. This is only about 0.07 
dex higher for Ay = 1. Errors on metallicities include 
the emission flux and Ay uncertainties only. The mean 
error, is 0.12 dex. Systematic uncertainties are generally 
higher than the measured errors for large metallicities. 
For 4 galaxies with an upper limit on H/3 or [OIII] flux, 
a lower or upper limit on the metallicity are provided, 
respectively. 

LCS used for the CFRS sample the calibrator pro- 
vided by McGaugh (1991) and Ay = 1. For consistency 
with GDDS, we have re-calculated metallicities, apply- 
ing Eq. [3 (Table 3). Errors in the metallicities include 
an arbitrary uncertainty on Ay of 0.3 mag. Also for 
this sample, galaxies with an upper limit on H/3 or [OIII] 
fluxes provide a lower or upper limit on the metallicity, 
respectively. Two galaxies with no line detection, have 
no information on the metallicity. Table 3 also lists old 
LCS metallicities for the upper branch. Our new mean 
metallicity and error are 12 + log(0/H) = 8.78±0.17 and 
0.07 dex, respectively (the LCS mean is 0.05 dex lower). 

Figure HHI shows 12 + log(0/H) as a function of red- 
shift, for the GDDS, CFRS and TKRS samples. The 
mean metallicity of all 261 galaxies is 12 + log(0/H) = 
8.72 ± 0.17, i.e. roughly solar 17 . Although the redshift 
range (0.3 < z < 1.0) corresponds to a time interval of 
4.3 Gyr, no redshift evolution is detected (i.e. no sig- 
nificant metallicity-redshift correlation). We emphasize 
that galaxies in this plot include all masses, therefore the 
presence of a mass-metallicity relation at these redshifts 
(see JFJ) may obscure any redshift evolution of metallicity. 

6. LUMINOSITY-METALLICITY RELATION AT Z ~ 0.7 

In Figure ITTI we show metallicity vs. B absolute mag- 
nitude (or B luminosity) for the GDDS (filled circles), 
CFRS (open squares) and the TKRS (crosses; KK04) 
galaxies. The GDDS+CFRS sample (79 galaxies) shows 
a 4<7 significant correlation. The linear bisector best fit 
(Feigelson & Babu 1992) is expressed by: 

12 + log(0/H) = (-0.280±0.045)Af B ab + (2.977±0.937) 

(6) 

The fit does not change much (although the significance 
is slightly lower) if Ay = 1 is assumed for the GDDS 
galaxies. 

In Figure ITTI we show this bisector fit (solid line) to- 
gether with the one found by KK04 for 177 TKRS galax- 
ies [12 + log(0/H) = -0.193M b ,ab + 4.900, long-dashed 
line]. The short-dashed line is the L-Z relation for the 
SDSS 53,000 z ~ 0.1 galaxies (Tremonti et al. 2004; here- 
after T04), which is very similar to the one found by 
Jansen et al. (2000) for the local Nearby Field Galaxy 
Survey (NFGS). 

KK04 note that galaxies in their sample have higher 
luminosities with respect to local galaxies with similar 
metallicities, suggesting a redshift evolution of the L-Z 
relation. This is confirmed by the GDDS+CFRS sam- 
ple, but our low-luminosity galaxies seem to have lower 
metallicities than the TKRS galaxies. This can only in 
part be due to the slightly different mean redshifts (0.71 

17 We assumed as a solar metallicity the value provided by Al- 
lende Prieto et al. (2001), which is 12 + log(O/H) = 8.69. 
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Fig. 11. — 12 + log(0/H) as a function of B absolute magnitude 
Mg Ag (or luminosity, upper horizontal axis) for the GDDS (filled 
circles), CFRS (open squares) and TKRS (crosses; KK04). The lin- 
ear bisector fits (see text) for GDDS+CFRS and TKRS are shown 
as solid and long-dashed lines, respectively. The short-dashed line 
is the luminosity-metallicity relation found for the SDSS sample at 
z ~ 0.1 (Tremonti et al. 2004). 
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Fig. 12. — Metallicity as a function of stellar mass for GDDS 
(filled circles) and CFRS galaxies (open squares). Error bars in- 
clude measurement errors but not possible systematic effects. The 
straight line is the linear bisector fit. The correlation is m ore than 
6cr significant. The metallicity scale is as in Figure ITTI and the 
mass scale spans three orders of m agni tudes, as in the luminosity 
scale of the upper a;-axis of Figure 1111 This way it is possible to 
visually compare the point scatter in the two figures. 



GDDS: Redshift Evolution of Mass-Metallicity Relation 



13 



for GDDS+CFRS and 0.62 for TKRS) because the dif- 
ference corresponds to only 0.5 Gyr. KK04 observe a 
steepening of the slope of the bisector fit from —0.16 to 
—0.24 (but this include the CFRS sample) going from 
z ~ 0.3 to z ~ 0.9. The difference might be related to 
the dust extinction correction, which can be a function of 
the stellar mass. While for the TKRS galaxies, EWs have 
been used which are almost insensitive to the dust cor- 
rection, for our emission fluxes we have assumed a con- 
stant extinction correction both for GDDS and CFRS. 
If instead the dust correction is higher for more massive 
(than generally brighter) galaxies, the slope of the bi- 
sector fit would be lower. We also note that the large 
uncertainties and the different selection criteria for the 
three samples can explain the discrepancy. 

7. MASS-METALLICITY RELATION AT Z ~ 0.7 

Stellar masses and metallicities have been derived for 
a total of 57 GDDS and CFRS galaxies. One galaxy, 
SA15-4272, was removed from the GDDS sample for the 
large uncertainty on the H/3 flux. The median and mean 
redshift of the remaining 56 galaxies is z — 0.69 and 0.71, 
for GDDS and CFRS respectively. The stellar mass and 
metallicity distribution (Figure I12f> are strongly corre- 
lated, at the level of 6.3cr. The linear bisector fit gives: 

12 + log(0/H) = (0.478±0.058)Af* + (4.062±0.579) (7) 

The significance of the correlation is almost the same 
(6.6<r) if the 4 galaxies in the turnover region (12 + 
log(0/-ff) < 8.4) are removed from the sample. The bi- 
sector slope becomes shallower, but consistent, with the 
slope for the whole sample (0.409 ± 0.050). 

The scatter around the fit is ~ 0.2 dex. The scatter 
for the luminosity- metallicity distribution of Figure ITT1 is 
65% higher. Results do not change much if Ay = 1 for 
the GDDS galaxies. The low-rcdshift half of the sample 
(z < 0.7) is not significantly differently distributed than 
the high-redshift half of the sample (z > 0.7). 

We compare this distribution with the M-Z relation re- 
cently found by T04 for the z = 0.1 SDSS galaxies. To do 
that, we need to convert their relation to one consistent 
with our choice of IMF and metallicity calibrator. 

T04 used the IMF provided by Kroupa (2001). The 
conversion from a Kroupa IMF to a BG03 IMF gives 
Af*,BG03 = A/*,koi/1-2, or a difference of 0.08 dex. 

T04 measured metallicities by modeling the emission 
lines and stellar continuum, using the approach described 
by Chariot et al. (in preparation), while for our sample 
we used the KK04 calibrator. A comparison between the 
two methods indicates that the T04 metallicities are at 
the most ~ 0.1 dex higher (this is mainly for the massive 
galaxies). The M-Z polynomial fit for the SDSS galaxies 
was recalculated using the KK04 calibrator (C. Tremonti, 
private communication) . 

The "converted" T04 M-Z relation (where the KK04 
metallicity and BG03 IMF are used) is: 

12+log(0/H) = -2.4412+2.1026 log A/*— 0.09649 log 2 M* 

The overall difference with the T04 polynomial is not 
very large. The new relation and the ±1<t dispersion are 
shown in Figure IT""1 



In the same figure, we show our z ~ 0.7 galaxies (as in 
Figure 11*2*1 but error bars are omitted). The majority of 
the galaxies are distributed below the z ~ 0.1 relation, 
and preferentially these are galaxies with M* < 10 10 M Q , 
suggesting that massive galaxies reach high metallicities 
much earlier than low massive galaxies. In other words, 
the mass-metallicity relation evolves from being steep at 
high redshifts to being flatter at low redshifts. This pic- 
ture of massive galaxies being more quiescent in recent 
times than less massive galaxies is also favored by recent 
results on the SFR density of the Universe calculated or 
measured for different galaxy stellar masses (Heavens et 
al. 2004; Juneau et al. 2005) according to which massive 
galaxies ceased forming stars (hence stopped enriching 
the ISM with metals) much earlier than low-mass galax- 
ies. Low-mass galaxies have a long-lasting star-formation 
activity, and significantly produce stars and metals even 
in recent times. Those galaxies in the mass-metallicity 
diagram of Figure lT""l will migrate with time from the left- 
bottom region to the middle region. A simple closed-box 
model which can explain this, is discussed in ii9.2l 

Our conclusion is not inconsistent with the relatively 
higher metallicities found for the generally massive LBGs 
(Shapley et al. 2004) and the star-forming galaxies at 
z ~ 2 of the K20 (de Mello et al. 2004). The rectan- 
gle in Figure 1131 represents the region occupied by the 
LBGs. Masses have been converted to match our IMF 
and include recycling into the ISM (if recycling is ig- 
nored, stellar masses are 1.4 times higher). Metallicities 
are also converted to match our metallicity calibrator. 
Shapley et al. (2004) used the N2 calibrator (Pettini & 
Pagel 2004), which gives, for those stellar masses, metal- 
licities ~ 1.4 to 2.3 times smaller, from the less massive 
to the most massive galaxy, than the KK04 calibrator 
(C. Tremonti, private communication). On average the 
difference is A [12 + log(0/H)] = -0.23. 

8. POSSIBLE SYSTEMATIC EFFECTS 

One of the puzzling results of the M-Z relation obtained 
for the star-forming galaxies by T04 and here, is that 
these galaxies appear to be generally metal rich for their 
stellar mass. For instance, the stellar mass of the Milky 
Way is ~ 5 x 10 10 M Q (Portinari et al. 2004). If its 
metallicity is approximately solar, 12 + log(0/H) ~ 8.7, 
this is 0.4 dex lower than that predicted by the M-Z 
relation of the SDSS (12 + log(0/H) - 9.1). 

First we have to consider that there are systematic 
differences between different metallicity calibrators (see 
M8.2I) . Moreover, it is well known that the metallicity 
decreases with the distance from the galaxy center (Gar- 
nett et al. 1997). According to Rolleston et al. (2000), 
in the Milky Way the oxygen abundance 18 goes from 
12 + log(0/H) = 9 to 8.6, at 6 kpc and 12 kpc from 
the center, respectively, with a radial gradient of —0.067 
dex kpc -1 . Although the quantitative result depends on 
the empirical method used to measure metallicity, and 
although some extragalactic studies indicate lower cen- 
tral metallicities (Kennicutt et al. 2003), it is likely that 
the metallicity is higher than the average if the slit aper- 
ture used to collect the galaxy signal is small. The fiber 
aperture for the SDSS sample (3") corresponds to a me- 

18 This is derived using absorption lines in the photosphere of 
B-type dwarf stars. 



14 



Savaglio, Glazebrook, Le Borgne et al. 



9.2 



9 - 



X 
o 



■2 8.6 



1.4 



5.2 - 



~i 1 1 r 




Fig. 13. — Metallicity as a function of stellar mass for GDDS 
(filled circ les) and CFRS (open squares). The points are as in 
Figure [l*2l but error bars are omitted. The green short-dashed line 
and hatched area indicate the polynomial fit and ±lcr dispersion, 
respectively, derived for the SDSS Z ~ 0.1 galaxies (T04). The 
blue rectangle is the region occupied by the LBGs at 2.1 < z < 2.4 
(Shapley et al. 2004) after correcting from the N2 to the R23 
calibrator. The black-solid and blue long-dashed lines are the SDSS 
polynomial shifted to the right to match the galaxy distributions 
at z ~ 0.7 and z ~ 2.3, respectively. 



dian diameter of 4.6 kpc (T04). At a distance of 4.6/2 
kpc from the center, according to Rolleston et al. (2000), 
the abundance of the Milky Way is 9.2, and higher if one 
takes the mean abundances in the inner region. Thus it is 
not surprising that SDSS galaxies have high metallicities. 

In the following subsections we show that other pos- 
sible systematic effects are not going to dramatically 
change our conclusions on the evolution of the M-Z rela- 
tion. 

8.1. Aperture effects 

If there is a gradient of the metallicity with the distance 
from the center (Vila-Costas & Edmunds 1992; Zaritsky, 
Kennicutt & Huchra 1994; Garnett et al. 1997; Rolleston 
et al. 2000), it is important to make sure that this effect is 
not dominating when comparing samples obtained using 
different slit apertures. 

We considered that the central part of a galaxy is more 
luminous than the outer part, and at high redshift is 
going to dominate the signal in the spectra. Kewley et al. 
(2005) have investigated the aperture effects as a function 
of redshifts, and found that the metallicity is generally 
not more than 0.1 dex larger than the average in the 
~ 4 kpc central region. This 4 kpc corresponds to an 
angular size < 0.8" for galaxies at z > 0.4. For this 
reason, the aperture effects are not strongly affecting the 
z ~ 0.7 results, because the region probed is generally 
large enough, and because residual offsets are affecting 
galaxies in a similar way. 

The SDSS fiber aperture corresponds to a median pro- 
jected diameter of 4.6 kpc (circular area is 16.6 kpc 2 ), 



with a range extending from 1 to 12 kpc 2 , from the low to 
the high redshift end. However, T04 estimated a modest 
(0.1 dex) metallicity change in the same redshift inter- 
val. The GDDS extraction aperture rectangle is 0.75 x 0.8 
arcsec 2 and corresponds to 3.9 x 4.2 kpc 2 and 5.9 x 6.3 
kpc 2 , from z = 0.38 to 0.96. The median for the sample 
is 5.6 x 5.9 kpc 2 . The dispersion of the physical area in 
the SDSS sample is higher than in GDDS, and the me- 
dian value is lower. It is hard to estimate precisely the 
mean difference between the two apertures due to the 
large aperture dispersion in the SDSS, we do not expect 
the aperture effect to be important. 

For the CFRS sample, the slit aperture is twice the 
GDDS aperture in the dispersion direction (1.3 arcsec, 
or 9.3 kpc at z = 0.7) and 20" in the perpendicular 
direction, so the galaxy area probed is larger than the 
SDSS area. However, following Kewley et al. (2005) and 
considering the small metallcity dispersion observed in 
the SDSS sample, we expect the effect to be small. 

Another related effect is the galaxy inclination. T04 es- 
timated 0.2 dex higher metallicity in fully face-on galax- 
ies than in fully edge-on galaxies, for fixed mass. For 
15 of our GDDS galaxies we obtained HST/ACS F814W 
images (Abraham et al., in preparation). One-third are 
face-on, 1/3 are edge-on and 1/3 are intermediate (Fig- 
ure EJ. Because this is a random effect, it would make 
the dispersion in the metallicity larger, but any signif- 
icant difference between the z ~ 0.1 and the z ~ 0.7 
samples would remain basically unchanged. 

8.2. Different metallicity calibrators 

A very important systematic effect to take into account 
is linked to the several metallicity calibrators used by 
different studies. Kennicutt et al. (2003) showed that 
the metallicity based on an electron-temperature sensi- 
tive method is on average a factor of 2.5 lower than that 
derived using the empirical method by Kewley & Dopita 
(2002). This factor becomes ~ 2, when compared with 
the calibrator used by us (KK04). 

We do not discuss here what the best method is, be- 
cause this is still controversial and is a strong function 
of what is observationally available (e.g., what set of 
emission lines are used). However, because of these dif- 
ferences, we have paid particular attention to compar- 
ing metallicities after converting to the same calibrator, 
namely, the KK04 one. This is the minimum requirement 
necessary to make the detection of the redshift evolution 
of the M-Z relation more robust. 

8.3. Dust extinction correction 

Another possible bias to consider is the dust extinction. 
As stated in i!3.2l dust extinction does not have a huge 
effect on the metallicity. This decreases by only 0.2 dex if 
Ay goes from to 3. We applied Ay = 2.1 and Ay = 1 
for GDDS and CFRS galaxies, respectively. The choice 
of using a lower extinction in the CFRS is justified by 
the higher luminosity of the CFRS galaxies, not entirely 
accounted for by the relatively higher stellar masses of 
the same sample with respect to the GDDS sample. In- 
deed, the CFRS is a R selected survey that used the 4m 
Canada-France-Hawaii Telescope, thus CFRS very likely 
could not identify highly extincted galaxies. The GDDS 
is a if selected survey that used a larger telescope, so 
suffers less from this bias. 
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Fig. 14. — The points show the displacement in mass applied to 
the SDSS M-Z relation (for which AlogAf = by definition) to 
match the GDDS+CFRS and LBG distributions. This describes 
how much more massive a high-z galaxy with a given metallic- 
ity would roughly be with respect to local galaxies with similar 
metallicities. The error bars indicate the ~ la mass disp ersi on at 
the median metallicity in the samples, as seen in Figure The 
small error bar for the LBGs is the effect of the small sample. The 
straight line is the linear correlation of the three points (Eq.l^t. 



Brinchman et al. (2004) found that massive (10 10 
Mq < M* < 10 11 Mq) star-forming galaxies on aver- 
age suffer 3 times (in magnitudes) more extinction than 
low mass galaxies (10 8 M Q < M* < 10 9 M ), i.e., oppo- 
site to what we have assumed. Even if there was no bias 
in the CFRS sample, assuming a larger extinction in the 
more massive galaxies would have a small effect on the 
correlation, but the scatter would be larger. 

9. MODELING THE MASS-METALLICITY RELATION 
EVOLUTION 

9.1. An empirical model 

We now attempt to derive an empirical model describ- 
ing the redshift evolution of the mass-metallicity relation. 
This will allow to predict the stellar mass (or the metal- 
licity) of a galaxy at a given redshift and metallicity (or 
stellar mass). 

Let us arbitrarily assume that at high redshift the 
shape of the SDSS relation is preserved, but it moves 
towards higher masses. We prefer to move the SDSS 
relation in the x direction instead of the y direction, be- 
cause this can better reproduce z ~ 0.7 GDDS and CFRS 
points in the M-Z plane. The best match is obtained if 
the M-Z relation is moved to the right by A log M ~ 0.47 
(black solid curve in Figure EJ- At higher redshift only 
the few galaxies studied by Shapley et al. (2004) have 
measured mass and metallicity, and these are represented 
by the rectangle in Figure El We have corrected metal- 
licities in order to take into account systematic effects. 
Shapley et al. (2004) have used the N2 calibrator, which 



saturates for high metallicities (Pettini & Pagel 2004). 
If we again assume that the shape of the M-Z rela- 
tion is preserved, the curve intercepts the locus of the 
LBG points for A log AT ~ 1.35. This mass displace- 
ment, the one at z ~ 0.7, and the SDSS point for which 
AlogM = 0, are shown in the plane tn - AlogM (tn 
is the Hubble time) of Figure ITU These three points are 
basically aligned along the straight line (also shown) 

A log Af (i H ) = -2. 0436 \ogt H + 2.2223 (9) 

in the log — log space. If we combine Eq. [5| with Eq. [H] 
we derive the general relation that gives the metallicity 
of a galaxy at a given stellar mass and Hubble time: 

12 + log(0/H) = -2.4412 + 2. 1026[logM* + 
-AlogM(iff)] - 0.09649[logM* - AlogM {t H )] 2 



(10) 



which is equivalent to: 



12 + log(0/H) = -7.5903 + 2.5315 log M„ 



-0.09649 log 2 M* 



5.1733 log t H - 0.3944 log 2 t H + 
-0.4030 logiff logM* 



fill 



Eq.^^is shown in Fienre HBI and predicts that low mass 
galaxies have a steeper redshift evolution in metallicity 
than bigger galaxies. Of course we have to consider the 
big uncertainty due to the way metallicities are mea- 
sured. However, Ea. II II describes in a simple and power- 
ful fashion the evolution of the metal cosmic abundance 
in the Universe for different mass bins. It predicts the 
mean metallicity of galaxies at a given mass and redshift, 
but does not give the time evolution of a single galaxy. 
This is done in the next section. 

9.2. A closed-box model 

A qualitative explanation of the redshift evolution of 
the M-Z relation (Figure^J is provided by a SFR which 
stops earlier in more massive galaxies than in low-mass 
galaxies (Heavens et al. 2004; Juneau et al. 2005). An- 
other way to describe this is a gas fraction (gas mass 
over total mass) available for star formation, that de- 
clines more rapidly in more massive galaxies, or a period 
of star formation that lasts longer in less massive galax- 
ies. 

Theoretical scenarios explaining these observations, 
must properly consider gas flows in both directions (in 
and out of the galaxy) . Galactic flows are mainly driven 
by the formation of massive stars which blow out gas 
from the galaxy, and by the gravitational potential of 
the galaxy (related to its mass), which on the one hand 
confines the gas in the galaxy, and on the other attracts 
metal-poor gas from the external IGM. The chemical evo- 
lution of galaxies over cosmic times was predicted by Pei 
& Fall (1995) using infall, outflow and closed-box (no 
flows) models. Larson (1974) described in his early work 
that low mass galaxies are metal poor due to the gas 
loss after supernova explosions. Ferrara & Tolstoy (2000) 
found later on in their models that indeed galaxies with 
gas mass M gas < 10 9 M lose mass in outflows. 
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Fig. 15. — 12 + log(0/H) as a function of stellar mass for constant redshifts (left panel), or as a function of Hubble time for constant 
stellar masses (right panel). These relations are obtained by extrapolating at all redshifts the z ~ 0.1 SDSS and our z ~ 0.7 M-Z relations. 
We emphasize that this empirical model is derived using the KK04 metallicity calibrator and BG03 IMF, and is different if other calibrators 
and IMFs are used. 



We tried to reproduce the M-Z relation evolution using 
the closed-box approach. Although the large dispersion 
in the M-Z point distribution can indicate that galactic 
flows play a role in the chemical state of galaxies, the 
closed-box model allows one to start to understand what 
basic assumptions are necessary to explain the observa- 
tions. We use the PEGASE modeling (Fioc & Rocca- 
Volmerange 1997, 1999) and assume that initially the 
galaxy is only made of gas with a given total mass M to t 
and zero metallicity. As the gas collapses and stars form, 
the total mass remains constant, and part of the gas is 
turned into stars. Both stellar mass and metallicity in- 
crease simultaneously as the galaxy evolves. The other 
key assumption is an exponential star formation SFR 
oc exp~'/ T , with an e- folding time r proportional to the 
total baryonic mass of the galaxy, according to the fol- 
lowing empirical relation: 

logr = alogM tot + 6 (12) 

For a < 0, the SFR declines more rapidly in more mas- 
sive galaxies. We tune this equation to reproduce the 
proposed M-Z relation evolution of Ea.lllland FigurefHfl 
This occurs for a = —0.88 and b = 9.42 (small insert in 
Figurc HBl) . The results of the model, assuming a redshift 
of formation zt = 3.5, are shown by the dotted lines in 
Figure OH they are not very sensitive to Zf, if this is in 
the range Zf = 3 — 1000. A galaxy with a stellar mass 
of 10 80 M at z = 2.3 will have a stellar mass of 10 88 
Mq at z = 0.1, i.e. a factor of 6 higher than 9.5 Gyr 
before (z — 2.3). Another galaxy with stellar mass 10 9 ' 7 
M© at z — 2.3, will only have doubled its stellar mass 
after the same time interval has passed. The increase in 
metallicity is a factor of 10 in the less massive galaxy, 
and only a factor of 3 in the more massive galaxy. This 
simple model, where r is a function of the mass, provides 



a good explanation of the observed results. Without any 
variation of r, it is not possible to account for the M- 
Z relation at z ~ 0.1, and all galaxies would have solar 
metallicity, regardless of mass and redshift of formation 
(if in the range zj = 3 — 1000). 

10. DISCUSSION 

The mass-metallicity (or the luminosity-metallicity) 
relation is a well known phenomenon identified in the 
local Universe in several surveys (Faber 1973; Lequeuex 
et al. 1979; Garnett & Shields 1987; Skilhnan, Kennicutt, 
& Hodge 1989; Zaritsky et al. 1994). Although any cor- 
relation is more difficult to detect at high redshift, due to 
the limited size of surveys, it is already clear that galax- 
ies with a given luminosity are more metal poor than 
similarly luminous galaxies in the local Universe (Kobul- 
nicky & Koo 2000; Pettini et al. 2001; Lilly et al. 2003; 
Liang et al. 2004; Shapley et al. 2004). Only recently it 
was possible to firmly establish that the relation between 
mass and luminosity persists at high redshift (z ~ 0.7; 
KK04), thanks primarily to the large sample of galax- 
ies available (about 200 from TKRS) and to the large 
interval in luminosity spanned. Such a relation at high 
redshift is displaced with respect to the local relation. 

All this is generally based on the B luminosity, a pa- 
rameter that poorly constrains the galaxy mass, because 
the B luminosity is easily extincted by dust, and repre- 
sents short-lived massive stars more than the bulk of the 
stellar mass. Stellar mass is a more meaningful galaxy 
parameter. 

The stellar mass and metallicity of 56 galaxies in our 
sample (0.4 < z < 1) are correlated (at ~ 6a signif- 
icance level). This relation is much stronger than the 
luminosity-metallicity relation for the same sample (for 
which the dispersion around the best fit is 65% higher), 
or for the TKRS sample at similar redshifts. 
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Fig. 16. — Closed-box model predictions of 12 + log(0/H) as 
a function of stellar mass (dotted lines), for redshift of formation 
Zf = 3.5, and for increasing values of the initial mass (from left to 
right). Results are not very sensitive to Zf, if this is in the range 
Zf = 3 — 1000. Open circles mark the three redshift epochs z = 2.3, 
0.7 and 0.1, from bottom to top. The dashed and solid lines are 
results from observation s at z = 2.3, 0.7 and 0.1 from bottom to 
top (as shown in Figure [T51 see text). These lines are solid in the 
mass range covered by observations. The small insert shows the 
c-folding time as a function of t he to tal galaxy mass adopted by 
the model, also described by Eg. 1121 for a = —0.95 and b = 10.6. 
The mass intervals marked by the diamonds are as in the dotted 
lines of the M-Z relation. 



This M-Z relation is displaced towards higher stellar 
masses, and/or lower metallicities, with respect to the 
z ~ 0.1 M-Z relation for SDSS star-forming galaxies. 
The rate at which mass and metallicity evolve with time 
can be estimated, also using results from the 7 z ~ 2.3 
LBGs studied by Shapley et al. (2004). Although the 
z ~ 2.3 sample is too small to establish any relation 
between mass and metallicity, we used it to normalize a 
hypothetical M-Z relation at z ~ 2.3. Such an additional 
information is very important, because it doubles the 
time interval: the Universe at z ~ 2.3 is 4.6 Gyr younger 
than at z ~ 0.7, which is 4.9 Gyr younger than at z ~ 
0.1. We assume that the shape of M-Z relation found 
at z ~ 0.1 does not change with time, and move it on 
the stellar-mass axis, to find the best match with the 
observed points at z = 0.7 and 2.3. 

We parametrize this displacement of mass and derive 
a general relation for the metallicity that is a function 
of the stellar mass and of the Hubble time (Eq. El and 
Figure IT5|l . With this generalized relation, we can sta- 
tistically predict the mass of a galaxy at a given redshift 
and metallicity. Ea. llll describes in a simple and powerful 
fashion the evolution of the cosmic metal abundance as a 
function of the galaxy stellar mass. It can also be used to 
test predictions of semi-analytic models (Somcrville, Pri- 
mack & Faber 2001), analytic models (Pei, Fall & Hauser 
1999), or numerical simulations (Nagamine et al. 2004). 

Any attempt to detect any cosmic chemical evolution 



would fail if this is done combining galaxies with very 
different stellar masses. Our empirical model indicates 
that massive galaxies evolve in terms of metallicity more 
slowly than less massive galaxies. The DLAs offer the 
opportunity to investigate consistently the chemical evo- 
lution for a particular class of galaxies over a large frac- 
tion of the cosmic time. It is now commonly accepted 
that DLAs do show chemical evolution (Prochaska et al. 
2003). This evolution can be compared with our predic- 
tions to derive the stellar masses of DLA galaxies. Of 
course, big uncertainties can lead to systematic errors. 
For instance, it is hard to estimate any systematic dif- 
ference between the metallicity derived from absorption 
lines and emission lines. Moreover, statistically, DLAs 
are tracing the outskirt of galaxies (Chen et al. 2005) 
where metallicity is generally lower than the metallicity 
in the central star-forming region (Ellison et al. 2005). 
Chen et al. (2005) estimated a metallicity gradient of 
0.041 ±0.012 dex kpc- 1 . 

We compared DLA observations to our empirical 
model, after applying to the model a constant negative 
correction in metallicity of 0.3 dex. Such an assumption 
is justified by assuming a typical impact parameter of 
the DLA from the galaxy center of 10 kpc, and consid- 
ering that our empirical relation is derived from galaxies 
observed over a region that is on average about 7 kpc 
across. We take the mean metallicity in redshift bins 
derived from a sample of 87 DLAs after dust depletion 
correction (Savaglio 2001). This is similar to the one 
more recently derived by Prochaska et al. (2003) using 
100 DLAs. The result is shown in Figure IT71 Remark- 
ably, the best fit predicts that the stellar mass of DLA 
galaxies is M* = 10 8 82±0 - 65 M , all the way from z ~ 0.5 
to 4.1. This suggests that typical DLAs are not originat- 
ing in very low-mass dwarf galaxies, but are intermedi- 
ate mass systems. The estimated stellar mass would be 

2 times higher (+0.3 dex) if the metallicity derived us- 
ing the i?23 parameter is overestimated by a factor of 2 
(Kennicutt et al. 2003; KK04). 

Figure ^| and Eq. ^] also suggests that local dwarf 
star-forming galaxies, with metallicities of ~ 1/10 solar 
(Aloisi et al. 2003; and references therein), have stellar 
masses of the order of 10 7 M©. In the distant Universe, 
the star- forming galaxies at z ~ 2.1 of the K20 with stel- 
lar masses (after correcting for BG03 IMF) of ~ 10 11 
M© (Fontana et al. 2004) would have a mean metallic- 
ity slightly over solar, in agreement with results derived 
using UV stellar features by de Mello et al. (2004). In 
a previous GDDS work (Savaglio et al. 2004) we have 
studied the ISM absorption in a sample of 13 galaxies at 
mean redshift z ~ 1.6 and found indication of high metal 
enrichment. The average stellar mass of this sample is 
M+ = 10 10 4 M©, which leads again, according to Eg. 1111 
to metallicity slightly above solar. 

For the LBGs studied by Pettini et al. (2001) at z = 

3 — 3.4, applying the KK04 calibrator and a modest dust 
extinction of Ay = 1, the reported emission line fluxes 
of 4 galaxies give 12 + log(0/H) = 8.4 - 8.8. At the 
redshift of the galaxies, Eq. II II gives a mean stellar mass 
of M* = 10 10 - 4 M©, i.e. about a factor of 2 less massive 
than the mean value in the z ~ 2.3 LBGs of Shapley et 
al. (2004). 

Massive LBGs are often used to derive the high-z SFR 
history of the Universe (the Madau plot; Madau et al. 
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Fig. 17. — Comparison between mean metallicity in DLAs for 
different redshift bins (Savaglio 2001), and our empirical model of 
Eq. [HI The best fit gives a stellar mass for DLAs of 10 8 ' 82 
(solid line) with a dispersion of 0.65 dex (dashed lines). In this 
figure we made a correction of —0.3 dex to the model to account 
for a metallicity gradient in DLAs (Chen et al. 2005). The best- 
fit mass would be 0.3 dex higher if the R23 metallicity, used to 
calibrate the empirical model, is underestimated by a factor of 2 
(Kennicutt et al. 2003). 



1996). The Madau plot is then used to derive the cosmic 
chemical evolution. A discrepancy with the DLA metal- 
licity evolution (which is lower and steeper) was noticed 
by Madau, Pozzetti, & Dickinson (1998). The discrep- 
ancy is likely because the metallicity evolution of dif- 
ferent galaxies with substantially different stellar masses 
are compared. This is confirmed by by Hopkins, Rao, 
& Turnshek (2005), who found that DLAs have a dom- 
inant role in the SFR density only in the late Universe 
(z < 0.6), but at higher redshift their contribution is 
much less important. We also note that even if the ana- 
lytic model of the cosmic metallicity evolution by Pei et 
al. (1999) resembles the DLA metallicity evolution, this 
alone does not necessarily prove that the cosmic chemical 
evolution is well represented by the DLA galaxies. 

Our empirical model can be tested using recent results 
on a DLA sample at 1.7 < z < 3.7 (Ledoux et al. 2005) 
showing a correlation between Fell absorption equiva- 
lent width (representing the galaxy dynamical mass) and 
metallicity. Metal rich DLAs have larger Fell EW than 
metal poor DLAs, a result which is reminiscent of the 
M-Z relation. If properly translated, their correlation at 
a mean redshift of z ~ 2.4 is consistent with our M-Z 
relation (C. Ledoux, private communication). 

The observed evolution of the M-Z relation can be re- 
produced by a closed-box scenario where the e-folding 
time is shorter for more massive galaxies. A star forma- 
tion more concentrated in time in massive galaxies than 
in less massive galaxies is consistent with the "downsiz- 
ing" scenario of galaxy formation. Other models cannot 



be ruled out by our simple approach. 

One shortcoming of our simple model is that it makes 
no attempt to consider the impact of galaxy-galaxy merg- 
ers on the evolution. This is known to evolve rapidly with 
redshift (Le Fevre et al. 2000). Mergers have the most 
dramatic effect when they are between galaxies of similar 
mass, and in many such cases the merger product will no 
longer be star-forming and hence not be included in these 
samples. The existence of a tight M-Z relation implies 
that an equal mass merger will shift galaxies to the right 
in Figure ^| the net effect will be to make the model 
tracks less steep. A model which included merging would 
then have to have a different form for r — M to t in order 
to compensate for this effect. A full treatment of merg- 
ing would require a detailed semi-analytic (e.g., Springel 
et al. 2005) or hydrodynamical model (Nagamine et al. 
2004) of galaxy formation which is beyond the scope of 
this work, though we believe the empirical relations we 
have discovered will greatly inform such a comparison. 
We note that it is much more straight-forward to com- 
pare stellar mass with metallicity as galaxies assemble in 
a model as one has to make less assumptions (e.g., about 
dust) than one does with a luminosity comparison. 

11. SUMMARY 

In this paper, we investigated for the first time the 
high-z mass-metallicity relation. We used galaxies from 
the GDDS and CFRS, for which stellar mass estimates 
were possible thanks to the available multi-band optical- 
NIR photometry (Crampton et al. 1995; Le Fevre et al. 
1995; Lilly et al. 1995; McCarthy et al. 2001; Chen et al. 
2002). 

The metallicity was measured using the R23 and O32 
over a sample of star- forming galaxies at 0.4 < z < 1. 
For the CFRS sample, we recalculated values measured 
by Lilly et al. (2003) using the recent calibrator provided 
by Kobulnicky & Kewley (2004) . Stellar mass and metal- 
licity are measured simultaneously for 56 galaxies. 

Thanks to the large mass range, M* = 10 82 — 10 108 
Mq, and the size of the sample, we unambiguously de- 
tect a 6cr significant correlation between stellar mass 
and metallicity at z ~ 0.7 (Eq. 0). Such a M-Z re- 
lation is displaced with respect to the M-Z relation at 
z ~ 0.1, towards higher stellar masses, or lower metal- 
licities. Galaxies at a given metallicity are more massive 
than galaxies in the local universe with similar metallic- 
ity. We attribute the origin of this to the redshift evolu- 
tion of the mass-metallicity relation. Although the shape 
of the M-Z relation depends to some extend on the choisc 
of the metallicity calibrator, our result on the evolution 
of the M-Z relation is robust, as all metallicities are es- 
timated using the same calibrator. 

The stellar mass is more tightly correlated to metallic- 
ity than the luminosity (the dispersion around the best 
fit is 1.6 times smaller). Moreover, the redshift evolu- 
tion of the M-Z relation is more significant than the red- 
shift evolution of the metallicity found in large samples 
of galaxies with different stellar mass. 

We find the M-Z evolution to be more rapid in lower- 
mass galaxies indicating that they are still actively being 
constructed. In contrast the more massive galaxies have 
already reached solar metallicity at z = 1 indicating that 
the bulk of their star-formation has completed. 

Based on the observed results, we derive an empirical 
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model for the evolution of the mass-metallicity relation 
as a function of redshift which works surprisingly well. 
From this, one can estimate the expected stellar mass of 
a galaxy at a given metallicity and redshift. If our model 
is correct, we predict that the stellar mass of DLAs is 
M* = 10 8 8±0 - 7 M at any redshift. According to the 
same model, the few LBGs at ~ 3.1 studied by Pettini et 
al. (2001) may have stellar masses of the order of 2 x 10 10 
Mq. The predicted masses are ~ 2 times higher if the 
T e -based metallicity calibrator is used (Kennicutt et al. 
2003). 

Our empirical model is nicely reproduced by a toy sce- 
nario of galaxy formation where the important prescrip- 
tion is a star formation that proceeds more rapidly in 
more massive galaxies, as described by the downsizing 
scenario for galaxy formation. Any more detailed model 
of galaxy formation must also reproduce the observed 
M-Z relation as a function of redshift. 

Future programs can test our predictions. For in- 
stance, by measuring metallicities of galaxies with known 
masses in the z ~ 1.5 Universe. Galaxies of the GDDS 
and K20 (Fontana et al. 2004; Glazebrook et al. 2004) 
would be a suitable testbed to this goal. To measure 
metallicities, NIR spectroscopy would be required, not 
an easy, but possible, task for those faint galaxies. Stellar 
masses of DLAs can directly be measured by obtaining 
K band photometry of QSO fields. This can more effi- 
ciently be done using the NIR capabilities of HST over 
\cw-z targets, to limit the confusion with the QSO PSF. 



Alternately ground-based Adaptive Optics may advance 
to the point where this is possible with larger telescopes. 
Finally, our proposed model of star formation history 
of galaxies, which is a function of mass, can be com- 
pared with the SFR history of the Universe in mass bins 
(Juneau et al. 2004; Heavens et al. 2004). 
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a Fluxcs, in 10 18 erg s 1 cm 2 , are not corrected for dust extinction. 

b Corrected for Balmer stellar absorption. 

c Corrected for dust extinction assuming Ay = 2.1. 

d Errors do not include systematic uncertainties. 

c Vega magnitudes from LCIRS. 

f Stellar masses in units of solar masses. 
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Table 2 

gdds composite emission line fluxes 



Line ratio 



Measured 



Theory a 



10 4 K 5000 K 



10 4 K 5000 K 



H7/H/3 

tW/H/3 

He/H/3 

H8/JL3 

[OII]/H/3 

[OIIL1A5007/H/3 

[OIIL1A4959/H/3 

[OIII]A4363/H/3 

[NcIII]/H/3 



0.292 ± 0.026 
0.180 ±0.017 
0.127 ±0.018 
0.035 ± 0.016 
2.280 ±0.081 
2.262 ± 0.088 
0.647 ± 0.042 

< 0.03 
0.208 ± 0.024 



0.470 0.458 3.47 ±0.64 3.28 ± 0.64 

0.262 0.250 1.91 ±0.48 1.67 ±0.48 

0.159 0.153 1.00 ±0.61 0.83 ± 0.61 

0.107 0.102 4.57 ± 1.89 4.37 ± 1.89 



Weighted mean 2.13 ± 0.32 1.92 ± 0.32 



a As expected from atomic physics, for two temperatures, and assuming no dust extinction. 
b Visual extinction, given the observed flux ratio and assuming two temperatures. 



GDDS: Rcdshift Evolution of Mass-Metallicity Relation 
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8.82 


± 0.10 


10.05 ± 0.40 


14 


0538 





810 












-21 


42 


8.76 


8.79 


± 0.03 




14 


0605 





837 












-20 


94 


8.78 


8.81 


± 0.05 




14 


0725 





580 


20 


39 


± 





23 


-19 


92 


9.00 


9.02 


± 0.03 


10.13 ± 0.25 


14 


0779 





580 












-20 


28 


8.98 


9.00 


± 0.05 




14 


0818 





901 












-22 


54 


8.86 


8.89 


± 0.08 




14 


0848 





664 


26 


07 


± 


9 


99 


-20 


37 


8.68 


8.73 


± 0.12 




14 


0972 





677 


20 


34 


± 





22 


-21 


56 


8.64 


8.69 


± 0.06 


9.92 ± 0.36 


14 


0985 





809 


21 


20 


± 





25 


-20 


76 


8.92 


8.94 


± 0.05 


9.97 ± 0.26 


14 


1087 





659 


21 


20 


± 





25 


-20 


63 


8.70 


8.72 


± 0.08 


9.53 ± 0.35 


14 


1126 





746 












-20 


71 


8.30 


8.53 


± 0.14 




14 


1189 





753 


21 


04 


± 





22 


-20 


86 


8.88 


> 


8.81 


9.84 ± 0.32 


14 


1190 





754 












-21 


99 


9.08 


> 


9.05 




14 


1258 





647 


21 


51 


± 





34 


-20 


31 


8.82 


8.87 


± 0.03 


9.59 ± 0.28 


14 


1386 





744 


19 


90 


± 





08 


-21 


67 


8.98 


9.01 


± 0.04 


10.21 ± 0.44 


14 


1466 





674 


23 


40 


± 


3 


70 


-20 


36 


8.30 


8.57 


± 0.12 





24 



Savaglio, Glazebrook, Le Borgne et al. 



Table 3 
- Continued 



ID 2 K M b ,ab* 12 + log(0/H) log M, 



(mag) (mag) LCS This papor b [M ] 



14 


9705 





609 












-21 


22 


8 


78 


8 


83 


± 





10 




22 


0274 





504 


19 


90 


± 





10 


-21 


60 


8 


52 


8 


53 


± 





11 


9.88 ± 0.39 


22 


0322 





915 


21 


12 


± 





31 


-21 


80 


8 


30 


8 


35 


± 





22 


9.93 ± 0.39 


22 


0417 


{) 


593 


21 


22 


± 





34 


-20 


30 


8 


80 


8 


86 


± 





08 


9.73 ± 0.35 


22 


0429 





624 


20 


49 


± 





17 


-20 


48 


8 


76 


8 


80 


± 





09 


10.07 ± 0.20 


22 


0576 





887 












-21 


24 


8 


70 


8 


72 


± 





04 




22 


0599 





886 












-21 


78 


8 


76 


8 


80 


± 





04 




22 


0770 





816 


22 


29 


± 





90 


-21 


47 


8 


30 


8 


60 


± 





12 




22 


0919 





472 












-20 


22 


8 


30 


8 


41 


± 





04 




22 


1119 





514 












-21 


87 


8 


86 


8 


90 


± 





08 




22 


1313 





817 


21 


16 


± 





32 


-21 


54 


8 


54 


8 


61 


± 





14 


9.92 ± 0.45 


22 


1350 





510 


21 


50 


± 





44 


-19 


88 


8 


80 


8 


84 


± 





06 


9.55 ± 0.40 


22 


1528 





665 


21 


84 


± 





60 


-20 


54 


8 


60 


8 


65 


± 





06 


9.39 ± 0.41 



a As reported in LCS. 

b Errors do not include systematic uncertainties. 



